clc;
clear;

load('AE_EDC_2.mat');

%提取各类条件
extent = EDC.extent;
traj = EDC.traj;
type = 1;

%初始化一个涡旋热量输送类
ETTEC = EddyThermalTransortClass('traj',traj,'extent',extent,'type',type);

%设置区域的间隔即将extent按照这个间隔进行子区域划分
ETTEC.intervalSet(1);

%首选需要将ETTEC类内的traj计算出trajThermal
ETTEC.trajThermalMeasure();

%计算热量输送
[MHT, ZHT] = ETTEC.eddyHeatTransportMeasure();

save('AE_ETTEC_2.mat','ETTEC','ZHT','MHT');